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(54) Design of a shielded gradient coil assembly for magnetic resonance imaging usingf inite 
element analysis 

(57) A method of designing shielded gradient coil 
assemblies (22) with a flared primary coil (60) for mag- 
netic resonance imaging systems is provided. The 
method includes generating a first continuous current 
distribution for a primary coil (60) using an inverse 
approach. The first continuous current distribution is 
confined within predetermined finite geometric bounda- 
ries of a surface defined by three dimensions and gen- 
erates a magnetic gradient field across an imaging 
region (14). The magnetic gradient field is constrained 
to predetermined values at specified spatial locations 
within the imaging region (14). The current distribution 
and magnetic field are converted into a stored energy 
and magnetic field domain where a finite element anal- 
ysis is performed to generate a second continuous cur- 
rent distribution for a shielding coil (62). The second 
continuous current distribution is confined within prede- 
termined finite geometric boundaries of a surface sur- 
rounding the primary coil (60). The second continuous 
current distribution generates a magnetic field which 
substantially cancels in an area outside a region defined 
by the shielding coil a fringe magnetic field generated by 
the first continuous current density. 
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[0001 ] The present invention relates to electro-magnetic shielding. It finds particular application in conjunction with 
gradient coils for magnetic resonance imaging apparatus and will be described with particular reference thereto How- 
ever, it is to be appreciated that the present invention also finds application in conjunction with systems which employ 
gradient magnetic fields and other applications in which electro-magnetic shielding is desired. 
[0002] In a magnetic resonance imaging system, gradient coil assemblies are commonly pulsed with electrical current 
pulses to produce magnetic gradients across the main magnetic field in the vicinity of an imaging region As an 
unwanted side effect, magnetic field gradients are produced which interact with external metallic structures such as the 
magnet cold sh.elds. the magnet dewar. and the like. The interaction generates eddy currents in the affected structures 
These eddy currents, in turn, generate eddy magnetic fields which have a deleterious effect on the temporal and spatial 
quality of the magnetic field in the vicinity of the imaging region and. hence, in the resultant image quality. 
[0003] The eddy current problem is often addressed by placing an active shielding coil between the primary gradient 
coil and the affected structure. The shielding coils are designed to substantially zero or cancel the magnetic field exter- 
nal to itself thereby preventing the formation of eddy currents in potentially vulnerable structures. 
[0004] Previously methods for production for magnetic gradients in magnetic resonance imaging systems consisted 
of winding discrete coils in a bunched or distributed fashion on an electrically insulating hollow cylindrical former and 
driving the coils with a current source of limited voltage. Conventional bunched coil designs include the Maxwell and the 
Modified Maxwell Pair for z-gradient production, and the Golay or Modified Golay (multi-arc) Saddle Coils for x and/or 
y-gradient production. Typically, these methods consisted of iteratively placing coil loops or arcs on the cylindrical 
former until the desired gradient strength, gradient uniformity, and inductance (related to stored energy) were achieved 
These previous designs were generally developed in a "forward approach" whereby a set of initial coil positions were 
defined (i.e.. the initial coil distribution), the fields and the inductance/energy calculated, and if not within particular 
design parameters, the coil positions would be shifted (statistically or otherwise) and results re-evaluated The iterative 
25 procedure continued until a suitable design was obtained. 

[0005] More recent methods of generating magnetic fields in magnetic resonance imaging systems utilize an "inverse 
approach." In the "inverse approach" method, the gradient magnetic field is forced to match predetermined values at 
specified spatial locations inside the imaging volume and a continuous current density is calculated which is capable of 
producing such a field. The "inverse approach" method assumes that the primary gradient coil has finite dimensions 
30 while those of the secondary or shield coil are left unrestricted (infinite). After the generation of continuous current dis- 
tributions for both the primary and the shield coils, an apodization algorithm is performed on the continuous current den- 
sity of the shield coil in order to restrain it to desirable dimensions. Following the modification of the shielding coil's 
continuous cunent. the Stream Function technique is employed in order to obtain discrete cunent patterns for both 
coils. Application of the Biot-Savart law to the discrete cunent pattern ensures that the discretization procedure was 
proper. This approach created generally more energy efficient gradient coil assemblies with higher gradient strengths 
and faster slew rates as compared to the "forward approach" method. 

[0006] In certain applications, for example imaging of the head, it is desirable to maintain the sweet spot of the gradi- 
ent coil assembly at a predetermined location corresponding to a particular anatomy of interest. However, due to a sub- 
ject's surrounding anatomy (the subject's shoulders in this case) the gradient coil assembly can not be appropriately 
centred. That is to say. when the gradient coil assembly's diameter is large enough to encompass a subject's shoulders 
it is too far radially from the subject's head, and conversely, when the diameter is sufficiently close to the head the 
shoulders interfere with an end of the coil such that the isocentre of the gradient coil is not properly aligned with the area 
of interest in the head. Shortening the length of the gradient coil so that the shoulders of a subject do not come into play 
reduces the usable imaging volume and deteriorates linearity. Attempts to adjust the sweet spot of the coil by means of 
45 a non-symmetric gradient coil assembly runs foul of high torques generated by the interaction of the main magnetic field 
with the currents of the gradient coil assembly and less desirable gradient fields. One approach used to address the 
issue, was to employ a symmetrical cylindrical gradient coil assembly with flared ends that could accommodate a sub- 
ject's shoulders. Added benefits of the flared design are increased patient access and reduced claustrophobia How- 
ever, heretofore shieldjng for this flared design has been lacking. 
so [0007] In one particular approach discussed in U.S. Patent No. 5.497.089 to Lampman et al.. an insertaWe cylindri- 
cal^ shaped gradient coil assembly is presented with flared ends. However, the gradient coil assembly is not shielded 
In this case, eddy cunent effects not suppressed by shielding prohibit the use of fast imaging sequences that are par- 
ticularly sensitive thereto. 

[0008] Another particular gradient coil assembly is presented in a paper by Schenck et al.. "Design Criteria for a 
55 Folded Gradient Coil". 5th ISMRM (Vancouver. Canada). April 1 997. pg. 1468. The paper presents a design methodol- 
ogy for a flared shielded gradient coil using the "forward approach". The design model presented assumed a series con- 
nection between the primary and shield coils such that they share the same current. This design resulted in inadequate 
shielding as compared to a traditional shielded design. 
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[0009] U.S. Patent No. 5,378,989 to Barber et al. presents yet another flared gradient coil assembly for use with open 
magnet systems. However, the flared portions are not at the ends of a cylindrical coil. Rather, the flared portions are 
near the isocentre. This prohibits the assembly from being used with a closed bore magnet Furthermore, the design 
restricts the flare to a 90° angle and is based on a "forward approach" which results in decreased performance. 

s [0010] For interventional procedures and like applications where patient access is desirable, it is advantageous to 
design the gradient shielding coil such that its dimensions do not exceed those of the flared primary gradient coil. In this 
manner, patient access can be maximized and the feeling of openness can reduce patient claustrophobia. However, in 
general, previous methods and prior art suffer the drawback that as the shielding coil length approaches that of the pri- 
mary coil, increased levels of eddy current effects within the imaging region deteriorate image quality. Conversely, when 

10 sufficient shielding is achieved, the dimensions of the shielding coil are substantially larger than those of the primary 
coil such that the level of patient access is encumbered and an increased level of patient claustrophobia is experienced. 
[001 1 ] In accordance with one aspect of the present invention, a method of designing shielded gradient coil assem- 
blies for magnetic resonance imaging systems is provided. The method includes calculating a first continuous current 
distribution for a three dimensional construction using an inverse approach. The first continuous current distribution is 

15 confined within predetermined finite geometric boundaries of the three dimensional construction and generates a mag* 
netic gradient field across an imaging region. The magnetic gradient field is constrained to predetermined values at 
specified spatial locations within the imaging region. Thereafter, a second continuous current distribution for a shielding 
coil is calculated using finite element analysis. The second continuous current distribution is confined within predeter- 
mined finite geometric boundaries of the shielding coil. The second continuous current distribution generates a mag- 

20 netic field which substantially cancels in an area outside a region defined by the shielding coil a fringe magnetic field 
generated by the first continuous current density. 

[0012] In accordance with another aspect of the present invention, a magnetic resonance scanner includes a main 
magnet for generating a main magnetic field through an examination region. The main magnet is arranged such that its 
geometry defines the examination region. A couch supports a subject to be examined within the examination region. A 

25 radio frequency coil is disposed adjacent the examination region for transmitting radio frequency signals into the exam- 
ination region and selectively exciting dipoles disposed therein. A radio frequency transmitter drives the radio frequency 
coil. A receiver receives magnetic resonance signals from resonating dipoles within the examination region. An image 
processor reconstructs an image representation from the received magnetic resonance signals for display on a human 
readable display. The scanner also includes a gradient coil assembly for generating substantially linear magnetic gradi- 

30 ents across the main magnetic field. The gradient coil assembly includes at least one cylindrical primary coil assembly 
with flared ends disposed about the examination region including an array of conductive coil loops arranged on a finite 
former such that a current density flowing thereon generates the substantially linear magnetic gradients. The gradient 
coil assembly also includes at (east one cylindrical secondary coil assembly disposed about the primary coil assembly 
between the primary coii assembly and the main magnet including an array of conductive coil loops arranged on a finite 

35 former such that a current density flowing thereon substantially cancels magnet flux density outside an area defined by 
the secondary coil assembly. The current density of the secondary coil assembly counters eddy currents that would be 
generated in the secondary coil assembly by the primary coil assembly. 

[0013] Ways of carrying out the invention will now be described in greater detail, by way of example, with reference 
to the accompanying drawings, in which: 

40 

FIGURE 1 is a diagrammatic illustration of a magnetic resonance imaging apparatus including a shielded gradient 
coil assembly designed in accordance with the present invention; 

FIGURES 2A and 2B are diagrammatic illustrations of a shielded gradient coil assembly designed in accordance 
with the present invention; 

45 FIGURE 3 is a three dimensional model of a portion of a shielded gradient coii assembly employed in the design 
thereof in accordance with the present invention; 

FIGURE 4 is a vector plot of the current distribution of the portion of the primary coil modelled in FIGURE 3 
designed in accordance with the present invention; 

FIGURE 5 is a vector plot of the current distribution of the portion of the shielding coil modelled in FIGURE 3 
so designed in accordance with the present invention; 

FIGURE 6 is a plot of the z component of the gradient magnetic field in Tesia versus distance along the x axis in 
meters in the region inside the primary coil at an axial position z = 0 m produce by the gradient coil assembly mod- 
elled in FIGURE 3 design in accordance with the present invention; 

FIGURE 7 is a plot of the z component of the fringe gradient magnetic field in Tesia versus distance along the z axis 
55 in meters at a radial position p = 0.485 m produce by the gradient coil assembly modelled in FIGURE 3 design in 
accordance with the present invention; and 

FIGURE 8 is a plot of the z component of the fringe gradient magnetic field in Tesia versus distance along the z axis 
in meters at a radial position p = 0.485 m produce by an unf lared shielded gradient coil assembly designed with a 
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truncated infinitely extended shielding current 

[0014] With reference to FIGURE 1 , a main magnetic field control 10 controls superconducting or resistive magnets 
1 2 such that a substantially uniform, temporally constant main magnetic field is created along a z axis through an exam- 
5 ination region 14. A couch (not illustrated) suspends a subject to be examined within the examination region 14. A mag- 
netic resonance echo means applies a series of radio frequencies (RF) and magnetic field gradient pulses to invert or 
excite magnetic spins, induce magnetic resonance, refocus magnetic resonance, manipulate magnetic resonance, spa- 
tially and otherwise encode the magnetic resonance, to saturate spins, and the like to generate magnetic resonance 
imaging and spectrography sequences. More specifically, gradient pulse amplifiers 20 apply current pulse to selected 
10 ones or pairs of gradient coil assemblies 22 to create magnetic field gradients along x, y, and z axes of the examination 
region 14. A digital radio frequency transmitter 24 transmit radio frequency pulses or pulse packets to a whole-body RF 
coil 26 to transmit RF pulses into the examination region. A typical radio frequency pulse is composed of a packet of 
immediately contiguous pulse segments of short duration which taken together with each other and any applied gradi- 
ents achieve a selected magnetic resonance manipulation. For whole-body applications, the resonance signals are 
15 commonly picked up by the whole-body RF coil 26. 

[001 5] For generating images of local regions of the subject, specialized radio frequency coils are placed contiguous 
to the selected region. For example, an inseiiable RF coil may be inserted surrounding a selected region at the isocen- 
tre of the bore. The inseiiable RF coil is used to excite magnetic resonance and receive magnetic resonance signals 
emitting, from the patient in the region being examined. Alternatively, the insertable RF coil can be used to only receive 
20 resonance signals introduced by whole-body coil RF transmissions. The resultant radio frequency signals are picked up 
by the whole-body RF coil 26, the insertable RF coil, or other specialized RF coils and demodulated by a receiver 30, 
preferably a digital receiver. 

[001 6] A sequence control circuit 40 controls the gradient pulse amplifiers 20 and the transmitter 24 to generate any 
of a plurality of multiple echo sequences such as echo planar imaging, echo volume imaging, gradient and spin echo 

25 imaging, fast spin echo imaging, and the like. For the selected sequence, the receiver 30 receives a plurality of data 
lines in rapid succession following each RF excitation pulse. Ultimately, the radio frequency signals received are 
demodulated and reconstructed into an image representation by a reconstruction processor 50 which applies a two 
dimensional Fourier transform or other appropriate reconstruction algorithm. The image may represent a planar slice 
through the patient, an array of parallel planar slices, a three dimensional volume, of the like. The image is then stored 

30 in an image memory 52 where it may be accessed by a display, such as a video monitor 54 which provides a human 
readable display of the resultant image. 

[0017] Heretofore, the conventional design for shielded gradient coil assemblies has generally been based upon the 
assumption that the total length of the shielding coil is unrestricted and thus its current density is assumed as infinitely 
extended. Consequently, as "infinite" shielding coils are not available for practical applications, the current density was 
35 subsequently apodized or truncated to fit on appropriate finite length elements. On the other hand, the technique 
employed herein involves the design of a shielded gradient coil configuration whereby the current for the shielding coil 
is initially restricted within the boundaries of a shielding coil having a finite length. 

[0018] With reference to FIGURES 2A and 2B. the geometrical configuration of the shielded gradient coil assembly 
22 designed such that the current density is initially contained within the finite boundaries of the shielding coil is shown. 

40 For this configuration, the length of the primary coil 60 is denoted as Lg, the axial length of the flared portion is denoted 
as L,. the length of the flared portion along its surface is defined as l^, and the length of the secondary coil 62 is denoted 
as Lb. The radius of the cylindrical portion of the primary coil 60 is denoted as a while the radius of the shielding coil 62 
is denoted as b. The flared portion of the primary coil is a conical surface with starting radius a and final radius a + e. 
Each combination of L^, Lf. and e defines a distinct conical surfaces. 

45 [0019] The design of a gradient coil where the z component of the magnetic field varies substantially linearly along 
the x direction (a transverse gradient coil or x-gradient coil) initially involves the design of the primary coil 60 based on 
an inverse approach methodology. Initially, the current of the primary coil is modelled to lie on a 2D cylindrical surface. 
For this x-gradient coil, the gradient magnetic field is anti-symmetric in the x direction around the geometric centre of 
the coil, while it is symmetric along the y and z directions. To generate such a gradient magnetic field, the analytical 

so expression of the electrical current for the primary coil can be written as: 

J*(0 = V\ *(<M)d^ + j z a (<f>.z)d z ]5(p-a) (1), 

where 5(p-a) is the restriction that confines the electrical current on a cylinder surface with radius a. The restriction to 
55 the primary coil length L^. the confinement of the current density on the cylindrical surface, the azimuthal and axial sym- 
metries for the jf and ;/ and the demand that the current density obeys the continuity equation provides the Fourier 
series expansion for both components of the current density around the geometric centre of the coil as follows: 



4 



EP 0 933 645 A2 



J + a (<M) = cos* £ for 



n=1 



(2). 



y z a (*,z) = sin* £ i^sin(/c n z) for |z|<s-^ 

n«1 n 



2 



(3). 
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where y^ a n are the Fourier coefficients, represents the total length of the primary coil, and 



2nn 



is since the electrical current cannot flow off the ends of the cylinder. Furthermore, both current components are zero for 

[0020] Expressing the z component of the magnetic field B 2 and the stored magnetic energy W m in terms of either 
one of the two components of the current density, the functional e in terms of W m and B z is constructed as: 



20 
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eU n a )-^m^ x j< B Arj)-B 2SC (rj)) 

7=1 



(4). 
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where are the Lagrange multipliers and B zS c represent the constraint values of the z component of the magnetic field 
at specified N points. Minimizing e, a quadratic function of the current, with respect to the current coefficients y^ a n 
obtains a matrix equation which j£ n satisfies: 



30 
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(5). 



where the evaluation of the Lagrange multipliers is done by way of the constraint equation. Inverting the matrix equa- 
35 tion, a solution for jf n is obtained and, hence, the solution for the current density is obtained as well. Replacing this 
expression back into the stored energy and magnetic field formulas, a final expression for the stored magnetic energy 
and the magnetic field is acquired in terms of the constraint points and the geometry of the system. In this manner, the 
form of the continuous current distribution along the azimuthal and axial directions tor the primary coil 60 is generated. 
[0021] In one embodiment for the design of the primary x-gradient coil, the radius a of the, cylinder is equal to 
40 0.342138 m and its total length l_a is restricted to 1.32700 m. Four Fourier coefficients were employed to specify the 
behaviour of the current density. Furthermore, three constraint points were employed to control the field behaviour 
inside a 45 cm imaging volume. Table 1 lists the constraint set used for designing trie primary gradient coil 60. 
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Pi 






1 


0.001 


0 


0.0000327 


2 


0.225 


0 


0.007766 


3 


0.001 


0.225 


0.0000262 



Table 1: Constraint set used for the design of the primary 
gradient coil. Values for p and z are in meters, values for 
B 2SC (2n) are in T. 



5 



EP 0 933 645 A2 



[0022] The first constraint point defines a gradient strength for the primary coil 60 of 32.7 mT/m, the second constraint 
point specifies a 5% linearity of the gradient field along the gradient x-axis up to a distance of 22.5 cm (the radial dis- 
tance of the 45 cm imaging volume) from the isocentre of the gradient field, and the third constraint point specifies a 
20% uniformity of the gradient field inside the 45 cm imaging volume. 

[0023] Through the use of these constraints and the application of the inverse approach methodology, the values are 
obtained for the Fourier coefficients for the current density, of the primary coil 60. which determines its continuous cur- 
rent distribution. Next, the continuous current distribution is discretized and the Biot-Savart law is applied to the discrete 
^current pattern in order to verify the validity of the earlier minimization. 

[0024] The continuous current distribution generated from the previous analytical method was then applied to the 
flared surface shown in FIGURES 2A and 2B. Having developed the primary coil 60, the next step is to shield the mag- 
netic field generated thereby with a finite length current distribution using a Finite Element Analysis. Due to symmetry 
conditions, only one portion needs to be modelled, with the appropriate boundary conditions that ensure correct current 
behaviour for the primary and shielding coils. With reference to FIGURE 3, a 3-D finite element model of the gradient 
coil assembly is illustrated. Again.due to symmetry, only one-haff the axial length, z from 0 to 



and one quarter of the angular section. <f> from 0° to 90°. is modelled. Additionally, in order to ensure the correct behav- 
iour of the magnetic field both inside and outside the gradient coil structures, the finite element model was extended 
axially to 5.0 m and radially to 2.5 m. At these locations, the vector potential was set to zero. The primary coil 60 is con- 
sidered as a layer with a cylindrical portion inner radius of 0.34096 m and an outer radius of 0.343320 m. This provides 
for a mean radius of 0.342138 m which corresponds to the radius a previously determined for the cylindrical portion of 
the primary coil 60. Hence, the primary coil 60 has a thickness equal to 0.00236 m and a half length, 



of 0.6635 m. The flared portion is defined by a cone with an axial length L, = 0.2m and a length along the surface of the 
cone 0.21 1 m. The mean starting radius of the cone is set to 0.342138 m to correspond with the radius a of the 
cylindrical section of the primary coil 60. The mean ending radius a + e = 0.410818 m and the thickness of the conduc- 
tor in the conical flared portion is set at 0.002364 m. 

[0025] The secondary coil has an inner radius of 0.432354 m and an outer radius of 0.433354 m while its half length. 



is 0.663504 m. In this embodiment, the half length of both the primary and secondary coils are chosen to be substan- 
tially identical. However, the finite element analysis code has been set in a parametric way which permits changes to 
the dimensions of the primary and secondary coils to fit within desired design parameters. Given the dimensions of the 
layer that corresponds to the primary coil 60. the current load using the Fourier coefficient expansion as it was previ- 
ously derived is applied thereto. FIGURE 4 illustrates the current pattern behaviour for the primary coil 60. To find the 
appropriate current pattern which will shield the primary coil's magnetic field and will flow within the boundaries of the 
secondary coil's layer, the finite element problem is approached as an eddy current problem. The material for the shield- 
ing coil's layer is taken to be copper with a resistivity of 1 .76x1 0" 8 O-m at 25° C. Performing a harmonic analysis with a 
frequency of 0.1 MHZ. or other appropriate high frequency, the solution for the current pattern of the shielding coil is 
obtained as shown in FIGURE 5. Note, the high frequency is chosen in order to avoid skin depth effects. In this manner, 
the shielding coil 62 is designed such that the fringe magnetic field generated by the primary coil 60 is substantially can- 
celled in a region outside the area defined by the secondary coil 62. Next the z component of the magnetic field along 
the gradient x direction is evaluated inside the region defined by the primary coil 60 to ensure the quality of the gradient 
magnetic field. FIGURE 6 plots the resulting total gradient (including the shielding coil's contribution) generated as the 
z component of the gradient magnetic field versus the distance along the x axis in the region inside the primary coil 60. 
The resulting total gradient strength is calculated to be 19.3 mT/m with an 8.8% linearity at a radial location of 0.2386 
m from the isocentre of the coil and a gradient strength of 21 mT/m at the isocentre. With reference to FIGURE 7. the 
shielding behaviour of the shielded gradient coil assembly 22 is evaluated at the region outside the shielding coil arid at 
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a radial location of p = 0.485 m. As FIGURE 7 illustrates, the net fringe magnetic field maximum occurs at the coil's iso- 
centre and has a value of approximately 12.4 ^T. 

[0026] With reference to FIGURE 8, an evaluation of the fringe magnetic field outside the shielding coil for a conven- 
tionally designed coil having a truncated infinite shield current as previously described and having no flared ends on the 

5 primary coil is illustrated. The shielding coil in this case is dimensionally similar to the one designed using the finite ele- 
ment analysis previously described but was designed with an infinitely extended current that was truncated to fit an axial 
half length of 0.8 m. At a radial location of p = 0.485 m, the coil designed with the truncated infinitely extending current 
has a fringe magnetic field maximum of 1 .05 mT which occurs at an axial distance of 0.68 m from the isocentre of the 
coil. In the case of a primary coil having flared ends, the shielding characteristics obtained from the truncation method 

10 are worse still. In comparison, the employment of the hybrid technique combining the inverse method and finite element 
analysis as described herein as opposed to the employment of the conventional truncation of an infinitely extended 
shielding current results in a reduction of the gradient coil's total fringe magnetic field by a factor of approximately 100. 
[0027] The design of the y-gradient shielded gradient coil assembly receives the identical treatment as that of the x- 
gradient coil merely rotated 90° about a central axis. Similarly, the z-gradient coil receives a like treatment with slight 

75 variation. Further, while the shielded gradient coil assembly herein has been described with reference to particular geo- 
metric boundaries and parameters, it is to be appreciated that various sizes, lengths, and geometric configurations of 
shielded gradient coil assemblies may be designed to meet the constraints of numerous applications. For example, the 
specified current patterns can be changed to produce either better linearity at the price of coil efficiency or greater effi- 
ciency at the price of linearity. The dimensions of the cylindrical gradient coil can be changed to be larger or smaller 

20 according to preferred applications. Another alternate embodiment includes the design of shielding coils where the 
length of the shielding coil is smaller than the length of the primary coil. As well, the gradient coil assembly may be 
designed such that net thrust forces generated as a result of the interaction of non-uniformities in the main magnetic 
field and currents flowing through the gradient coil assembly are balanced for each primary and shielding coil set. Fur- 
thermore, it is to be appreciated that while described with reference to magnetic resonance imaging systems having 

25 central bore examination regions, the design technique is also applicable to open type magnetic resonance imaging 
apparatus having axial and/or vertically directed fields. 

[0028] One advantage of the present invention resides in the ability to design shielded gradient coil assemblies for 
short magnet geometries which improve patient access and reduce patient claustrophobia. Another advantage of the 
present invention resides in providing electro-magnetic shielding for a three dimensional current distribution. Another 
30 advantage of the present invention is reduced eddy current effects in the vicinity of the imaging region. Another advan- 
tage of the present invention resides in the production substantially linear magnetic gradients with relatively high 
strengths and increased slew rates. Another advantage of the present invention is that it permits the calculation of cur- 
rent patterns for a finite shielding coil without the use of current truncation or apodization mechanisms. 

35 Claims 

1 . A method of designing shielded gradient coil assemblies for magnetic resonance imaging systems comprising: (a) 
calculating a first continuous current distribution for a three dimensional construction using an inverse approach 
such that the first continuous current distribution is confined within predetermined finite geometric boundaries of 

40 the three dimensional construction and generates a magnetic gradient field across an imaging region, the magnetic 
gradient field constrained to predetermined values at specified spatial locations within the imaging region; (b) cal- 
culating a second continuous current distribution for a shielding coil using finite element analysis such that the sec- 
ond continuous current distribution is confined within predetermined finite geometric boundaries of the shielding 
coil, such that the second continuous current distribution generates a magnetic field which substantially cancels in 

45 an area outside a region defined by the shielding coil a fringe magnetic field generated by the first continuous cur- 
rent density. 

2. A method as claimed in claim 1 , further including: discretizing the first and second continuous current distributions 
to generate first and second winding patterns; applying the first winding pattern to a three dimensional primary coil; 

so and applying the second winding pattern to the shielding coil. 

3. A method as claimed in claim 1 or claim 2 wherein: the first winding pattern lies on a cylindrical surface with out- 
wardly flared ends; and the second winding pattern lies on a cylindrical surface. 

55 4. A method as claimed in claim 3, wherein the shielding coil has an axial length substantially equal to or less than 
that of the primary coil. 

5. A method as claimed in any one of claims 1 to 4, wherein the shielding coil is made of a conductive material and a 



7 



EP0 933 645 A2 

harmonic analysis at high frequency is performed to obtain the second continuous current distribution. 

6. A gradient coil assembly for use in a magnetic resonance scanner comprising: a main magnet (12) for generating 
a main magnetic field through an examination region (14), the main magnet (1 2) being arranged such that its geom- 
etry defines the examination region (14); a couch which supports a subject to be examined within the examination 
region (14); a radio frequency coil (26) disposed adjacent the examination region (14) for transmitting radio fre- 
quency signals into the examination region (14) and selectively exciting dipoles disposed therein; a radio frequency 
transmitter (24) for driving the radio frequency coil (26); a receiver (30) which receives magnetic resonance signals 
from resonating dipoles within the examination region (14); an image processor (50) which reconstructs an image 
representation from the received magnetic resonance signals for display on a human readable display (54); the gra- 
dient coil assembly (22) being suitable for generating substantially linear magnetic gradients across the main mag- 
netic field, and the gradient coil assembly including; at least one cylindrical primary coil assembly (60) with flared 
ends disposed about the examination region (14) including an array of conductive coil loops arranged on a finite 
former such that a current density flowing thereon generates the substantially linear magnetic gradients- and at 
least one cylindrical secondary coil assembly (62) disposed about the primary coil assembly (60) between the pri- 
mary coil assembly (60) and the main magnet (12) including an array of conductive coil loops arranged on a finite 
former such that a current density flowing thereon substantially cancels magnet flux density outside an area 
defined by the secondary coil assembly (62) 

7. A gradient coil assembly as claimed in claim 6, wherein the formers are hollow cylindrical tubes, the former for the 
primary coil assembly (60) having flared ends, arranged such that the examination region (14) is defined inside the 
former of the primary coil assembly (60) with the former of the primary coil assembly (60) positioned inside the 
former of the secondary coil assembly (62) and an axial length of the secondary coil assembly (62) is substantially 
equal to or less than that of the primary coil assembly (60). 

8. A gradient coil assembly as claimed in claim 6 or claim 7, wherein a maximum magnetic flux density outside the 
area defined by the secondary coil assembly (62) is less than 50 jiT. 

9. A gradient coil assembly as claimed in any one of daims 6 to 8. wherein the current density at axial ends of the 
cylindrical secondary coil assembly (62) is great©- that at a centre. 

10. A gradient coil assembly as claimed in any one of claims 6 to 9, wherein the gradient coil assembly (22) includes 
three primary coil assembly (60) and secondary coil assembly (62) pairs for generating substantially linear mag- 
netic gradients along three mutually orthogonal axes. 

1 1 . A gradient coil assembly as claimed in any one of claims 6 to 1 0. wherein the current density of the secondary coil 
assembly (62) is arranged to counter eddy currents that would be generated in the secondary coil assembly (62) 
by the primary coil assembly (60). 

12. A magnetic resonance scanner having a gradient coil assembly as claimed in any one of claims 6 to 11. 
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(54) Design of a shielded gradient coil assembly for magnetic resonance imaging usingfinite 
element analysis 

(57) A method of designing shielded gradient coil 
assemblies (22) with a flared primary coil (60) for mag- 
netic resonance imaging systems is provided. The 
method includes generating a first continuous current 
distribution for a primary coil (60) using an inverse 
approach. The first continuous current distribution is 
confined within predetermined finite geometric bounda- 
ries of a surface defined by three dimensions and gen- 
erates a magnetic gradient field across an imaging 
region (14). The magnetic gradient field is constrained 
to predetermined values at specified spatial locations 
within the imaging region (14). The current distribution 
and magnetic field are converted into a stored energy 
and magnetic field domain where a finite element anal- 
ysis is performed to generate a second continuous cur- 
rent distribution for a shielding coil (62). The second 
continuous current distribution is confined within prede- 
termined finite geometric boundaries of a surface sur- 
rounding the primary coil (60). The second continuous 
current distribution generates a magnetic field which 
substantially cancels in an area outside a region defined 
by the shielding coil a fringe magnetic field generated by 
the first continuous current density. 
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